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Abstract 

The adiabatic self-consistent collective coordinate (ASCC) method is applied to 
the pairing-plus-quadrupole (P + Q) model Hamiltonian including the quadrupole 
pairing, and the oblate-prolate shape coexistence phenomena in proton-rich nuclei, 
68 Se and 72 Kr, are investigated. It is shown that the collective path connecting 
the oblate and prolate local minima runs along a triaxial valley in the {[3, 7) plane. 
Quantum collective Hamiltonian is constructed and low-lying energy spectra and 
E2 transition probabilities are calculated for the first time using the ASCC method. 
Basic properties of the shape coexistence/mixing are well reproduced. We also 
clarify the effects of the time-odd pair field on the collective mass (inertial function) 
for the large-amplitude vibration and on the rotational moments of inertia about 
three principal axes. 

1 Introduction 

In proton rich nuclei around 68 Se and 72 Kr, oblate-prolate shape coexistence phenomena 
are observed [1—4]. From the mean-field point of view, the oblate and prolate shapes 
are mixed by the many-body tunneling through the potential barrier between the two 
local minima in the potential energy landscape. To describe the shape mixing dynamics, 
we have to determine the collective degrees of freedom in which direction the large- 
amplitude shape dynamics takes place. The adiabatic self-consistent collective coordi- 
nate (ASCC) method [5] has been proposed as a microscopic theory of large-amplitude 
collective motion such as the shape coexistence phenomena. This theory is based on the 
time-dependent Hartree-Fock-Bogoliubov (TDHFB) theory and enables us to extract 
the collective degrees of freedom in a self-consistent way. In the previous work [6], we 
have solved the pairing-plus-quadrupole (P + Q) Hamiltonian by means of the ASCC 
method for 68 Se and 72 Kr nuclei, and successfully extracted the one-dimensional collec- 
tive path connecting the oblate and prolate local minima. Since the two local minima 
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are mainly connected by the triaxial degrees of freedom, we have suggested that the 
triaxial deformation plays an essential role in the shape coexistence dynamics of these 
nuclei. 

In this presentation, we report results of the first calculation of low-lying energy 
spectra and E2 transition probabilities by means of the ASCC method. We derive 
the quantum collective Hamiltonian that describes the coupled collective motion of the 
large-amplitude vibration responsible for the oblate-prolate shape mixing and the three- 
dimensional rotation of the triaxial shape. The calculation yields the excited prolate 
rotational band as well as the oblate ground-state band. It also indicates that the shape 
mixing decreases as the angular momentum increases. 

We also clarify the effect of the time-odd mean-field on the collective mass. This effect 
is ignored in the cranking mass [7], but is included in the ASCC mass. The time-odd 
mean-field effect generated by the particle-hole residual interaction was investigated in a 
few decades ago [8], but those generated by the pairing interaction has not been discussed 
so far. Quite recently we have shown, using the schematic model Hamiltonian [9], that 
the time-odd component associated with the quadrupole-type pair field significantly 
increases the collective mass. In the present calculation, we thus include the quadrupole 
pairing interaction to the P + Q Hamiltonian to clarify the effect of the time-odd pair 
field on the collective mass and rotational moments of inertia. 



We first recapitulate the basic equations of the ASCC method. The moving-frame TD- 
HFB state \4>{q,p)) = e %p< ^^ \<j)(q)) is written in terms of the collective coordinate q, 
and collective momentum p. The TDHFB state is expanded in terms of the collective 
momentum p under the adiabatic assumption. The collective Hamiltonian is expanded 
up to the lowest order in p as 



where the V(q) and B(q) represents the collective potential and the inverse collective 
mass. We add the rotational energy term with three moments of inertia Ji(q) to the 
collective Hamiltonian in order to take into account the rotational motion. 

The basic equations of the ASCC method are derived from the adiabatic expansion 
of the time-dependent variational principle, and are summarized as 



2 The ASCC method 




(i) 



8(<t>{q)\H M {q)\<Kq))=U, 



(2) 



5 (0(g)! [H M (q), Q(q)] ~ -B(q)P(q) \<j>(q)) = 0, 



(3) 
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5{<t>(q)\[H M (q),~P(q))-C(q)Q(q) 

Equation (J2]) is called the moving-frame HFB equation, while Eqs. ([3j) and © are 
called the moving-frame QRPA equations, where Q(q) and P(q) denote the infinitesimal 
generators of the collective path, 

H M (q) =H-J2 \ {T) (q)N^ - ^f(q)Q(q), (5) 

T 

is the Hamiltonian in the moving-frame, and C(q) = + 2B(q) 7)f f>q" * s ^ e l° ca l 
stiffness parameter. The operators H, and the quantity \^(q) represents the 

microscopic Hamiltonian, the particle number operators, and the chemical potential, 
respectively. The QRPA eigenmodes about the HFB equilibrium points, as well as static 
solutions of the HFB equation, are always special solutions of the ASCC equations. 
Therefore, the collective path can be constructed by means of the local progression 
procedure starting from a static HFB state. The collective path at q + 5q is calculated 
using the constraint {(p(q + 5q) \ Q(q) \4>(q + 5q)) = 5q derived from the canonical variable 
conditions, if the collective path at q is already known. Repeating this procedure, we 
can find the collective path and obtain all quantities appearing in the ASCC equations 
and collective Hamiltonian except the rotational moments of inertia. Three rotational 
moments of inertia are evaluated by extending the QRPA equations for rotation at the 
HFB minima to the general HFB states \(f)(q)) on the collective path, 

5(<t>(q)\ [H M (q),i%(q)\ - jr\q)J % \<f>(q)) = 0, (6) 

where ^?i(q) and Jj are the angle and angular momentum operators in the principal 
frame. 

In the present calculation, we use the P + Q + quadrupole pairing Hamiltonian as 
a microscopic Hamiltonian H. We adopt the same single-particle energies and the P + 
Q interaction strengths as in Ref. 6 and the self-consistent quadrupole-pairing strength 
G| clf proposed by Sakamoto and Kishimoto [10]. 

3 Collective Path 

For both 68 Se and 72 Kr, we have found that the lowest HFB state possesses the oblate 
shape, while second lowest HFB state the prolate shape. We start from the oblate state 
{q = 0) and determine the collective path connecting the two local minima. Figure [T] 
shows the collective path projected onto the (/3,7) potential energy surfaces. The path 
connects the two local minima via the triaxially deformed region. Figure [2] shows various 
quantities appearing in the collective Hamiltonian and the basic equations of the ASCC 
method. We define the scale of collective coordinate q by setting B(q)^ 1 = 1. The 
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collective mass with respect to the conventional (/?, 7)- coordinates is then given by 

M(q) = ^(d(5/dq) 2 + f3 2 (dj / dq) 2 . We have found that the collective mass M(q) 
and three rotational moments of inertia Ji{q) are enhanced by the time-odd pair field 
generated by the quadrupole pairing. 



4 Energy Spectrum and Transition Probabilities 

We requantize the collective Hamiltonian ([T]) and solve the collective Schrodinger equa- 
tion 



(7) 



to obtain the excitation energies and collective wavefunctions. The collective wave func- 
tions iM,k(q,&) takes the following form: 



*IM,k{q, = E ®lK,k{q) < n\IMK >, (8) 
K=0 

where &iK,k{q) an d < Sl\IMK > denote the vibrational wavefunctions and the rotation 
matrices, respectively. The conventional boundary conditions and symmetry require- 
ments for solving the Bohr collective Hamiltonian [11] are adopted. 

Figures display the energy spectrum and B(E2) values for 68 Se and 72 Kr. The 
calculation yields the excited prolate rotational band as well as the oblate ground-state 
band. It is seen that the inter-band E2 transitions are weaker than the intra-band 
E2 transitions, indicating that the oblate-prolate shape coexistence picture holds. The 
result of calculation also indicates that the oblate-prolate shape mixing decreases as 
the angular momentum increases. Thus, the basic patterns of oblate-prolate shape 
coexistence are successfully reproduced in the calculation. The calculation suggest the 
existence of excited + state in 68 Se, which is not yet found in experiment. By comparing 
two calculations, (G2 = 0) and (G2 = C| elf ), we found that the energy of the 
state is much lowered by including the quadrupole pairing interaction. This is due 
to the enhancement of the collective mass. The energies of the other members of the 
rotational bands are also lowered, because the quadrupole pairing also enhances the 
rotational moments of inertia. Experimental data for B(E2) are not available except 
5(E2;2+ -> 0^)=1000 e 2 fm 4 in 72 Kr [4]. Theoretical B(E2) values are calculated 
using a polarization charge e po i = 0.881 that reproduces this data, so that only relative 
magnitudes should be regarded as theoretical estimates. 

In order to estimate the oblate-prolate shape mixing in quantum eigenstates, we 
define the oblate and prolate probabilities as follows: 

P oh {I,k)= dqJ2\®lK,k(<l)\ 2 i Ppro(I,k)= dqJ2\®lK,k(l)\ 2 , ( 9 ) 

•^<7min K=0 1° K=0 
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where we assume that g m j n < q Q ^ < qo < q pTO < q m &x- The "boundary" between the 
oblate and the prolate regions is set to the top of the potential barrier or at 7 = 30°. 
Figure U] shows these probabilities for 68 Se and 72 Kr. It is clearly seen that the shape 
mixing rapidly decreases as the angular momentum increases. 

5 Summary 

For the first time, we have reported the result of calculation of low-lying energy spectra 
and E2 transition probabilities in 68 Se and 72 Kr by means of the ASCC method. We have 
derived the quantum collective Hamiltonian that describes the coupled collective motion 
of the large-amplitude vibration responsible for the oblate-prolate shape mixing and 
the three-dimensional rotation of the triaxial shape. The calculation yields the excited 
prolate rotational band as well as the oblate ground-state band. It also indicates that 
the oblate-prolate shape mixing decreases as the angular momentum increases. 

It is surprising that basic pattern of the shape coexistence/mixing phenomena is well 
reproduced using the one-dimensional collective path running on the two-dimensional 

7) plane. Speaking more exactly, this collective path is self-consistently extracted 
from the huge dimensional TDHB manifold. Namely, the result of calculation indicates 
that the TDHB collective dynamics of the shape coexistence phenomena in these nuclei 
is essentially controlled by the single collective variable microscopically derived by means 
of the ASCC method. 

We have also shown that the time-odd pair field enhances the collective mass of the 
large-amplitude vibrational motion and the rotational moments of inertia. This finding 
is important in understanding the shape coexistence dynamics, because, together with 
the collective potential energy, these inertial functions associated with the collective 
kinetic energies determine the degree of localization of the collective wave function in 
the (/?, 7) plane. 
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Figure 1: Collective paths for 68 Se (upper) and 72 Kr (lower) calculated using the P + 
Q Hamiltonian including the quadrupole pairing interaction. The paths are projected 
onto the (/?, 7) potential energy surface. The dots in the figures indicate the HB local 
minima. The equipotential contour lines are drawn every 100 keV. 
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Figure 2: The monopole pairing gaps, Aq (q), the quadrupole pairing gaps, A^ (q) 
and A^iq), the collective potential V(q), the collective mass M(q), the rotational 
moments of inertia Ji(q), the lowest two moving-frame QRPA frequencies squared, 
^(q) = B(q)C(q), and the quadrupole deformations, (3(q) and 7(9), are plotted as 
functions of 7(g) for 68 Se (upper) and q for 72 Kr (lower). Results of two calculations 
using the P + Q Hamiltonian with (G2 = G| elf ) and without (G2 = 0) the quadrupole 
pairing interaction are compared. 
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Figure 3: Excitation energies and B(E2) values for low-lying states of 68 Se and 72 Kr 
calculated by the ASCC method. In the left (middle) panel, the quadrupole pairing is 
ignored (included). Experimental data [1-4] are displayed in the right panel. The B(E2) 
values are given in parentheses beside the arrows in unit of e 2 fm 4 . 
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Figure 4: The oblate and prolate probabilities evaluated for individual eigenstates in 
68 Se (left) and 72 Kr (right). For each nucleus, the left (right) panel shows values for 
the lowest (the second lowest) state of each angular momentum. The open (closed) 
symbols indicate the oblate (prolate) probabilities. The probabilities defined by setting 
the boundary at the barrier top (7 = 30°) are shown by squares (circles) In each figure, 
the upper (lower) panel shows the probabilities calculated without (with) including the 
quadrupole pairing interaction. 



10 



